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Abstract. In this paper we are investigating the long time behaviour of the solution of a mutation 
competition model of Lotka-Volterra's type. Our main motivation comes from the analysis of the 
Lotka-Volterra's competition system with mutation which simulates the demo-genetic dynamics of 
diverse virus in their host : 



dt 



K 



■^f^ijiVj -Vi). 



In a first part we analyse the case where the competition terms 'J'i are independent of the vims type 
i. In this situation and under some rather general assumptions on the functions the coefficients 
Ti and the mutation matrix fiij we prove the existence of a unique positive globally stable stationary 
solution i.e. the solution attracts the trajectory initiated from any nonnegative initial datum. More- 
over the unique steady state v is strictly positive in the sense that Vi > for all i. These results 
are in sharp contrast with the behaviour of Lotka-Volterra without mutation term where it is known 
that multiple non negative stationary solutions exist and an exclusion principle occurs (i.e For all 
i ^ io,Vi =0 and Vi^ > ). Then we explore a typical example that has been proposed to explain 
some experimental data. For such particular models we characterise the speed of convergence to 
the equilibrium. In a second part, under some additional assumption, we prove the existence of a 
positive steady state for the full system and we analyse the long term dynamics. The proofs mainly 
rely on the construction of a relative entropy which plays the role of a Lyaprmov functional. 
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1. Introduction 

In this paper we are investigating the long time behaviour of the solution of some models that 
have been recently used in epidemiology. Our analysis focuses on a Lotka-Volterra competition 
system with mutation which basically simulates the demo-genetic dynamic of a genetically di- 
verse virus population in its hosts, highlitrng the numerous links existing between ecological and 
withrn-host infection dyanmics. Such type of model has been proposed to explain some exper- 
imental data e.g. llT6ll22l |37|| . To be more specific the demo-genetic dynamic is modelled by N 
ordinary differential equations which simulates at host scale the dynamics of Vi (t) the number of 
virus particles of genotype i at time t : 



(1.1) — ^ 
^ ' dt 



N 



where r, , K and /iy represent respectively the growth rate for each genot3'pe, the total carrying 
capacity of the host and a nonnegative matrix characterising the rate of mutation from a virus of 
genotype i to a virus of genotype j. For each i, ^'^(u) : ^ M is a locally Lipschitz application 
describing the intensity of the interaction between a virus of genotype i with all its competitors. 

The mutation term of the system can also be interpreted as a dispersal term. Indeed, mutation 
naturally corresponds to dispersal into the discrete space of genotype. But the mutation term can 
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also handle dispersal between physical patches. With this in mind, the above system of equation 
can also be used to model the demo-genetic dynamic of a diverse virus population in structured 
hosts, each host tissue types being virus "habitats" connected to each others by dispersal via fluid 
flow (e.g. ||38l ). Ways of derive results for this interpretation are discussed in the biological 
comments subsection. 

In what follows we will always make the following assumptions on Ti, and /i^ 

' For all i, ri > 0, 

The matrix (yUij ) is nonnegative symmetric and irreducible 

eC,ti(M^,M),vI/,(0) = 

is monotone increasing with respect to the natural order of 

Furthermore we will assume that for all i there exists positive constants Ri,ki, Ci with fc^ > 
so that the function vj/j satisfies for all v e + \ Qr^ (0), 

ki 



(1.2) 



N 



(1.3) 



where Qr. (0) denotes the ball of radius Ri and centred at associated to the norm. 
A typical example of model satisfying our assumption is given by 



dvijt) 
dt 



I 4 



Mv 



with 



M 



m(i - m) 
m(i - m) 



m(i - m) (1 - - 

where is a parameter giving the point mutation rate per replication cycle and per nucleotide. 
This mutation matrix corresponds to a viral population composed of 4 variants differing only 
by one or two substitutions involved in adaptative proprieties (e.g. pathogenicity). The interac- 
tions terms ^'i :— Pij""] handle a wide range of possible inter-specific (rnter-variants) 
competition rates between any pairs of virus variants. 

This particular structure has been used recently to model the adaptation of plant virus to re- 
sistance genes ||2T1 122| . This particular form of competition is commonly used to model virus 
evolution, lll6ll32ll37ll . Without the mutation's matrix (/lij), the system l|l.H reduces to a classical 
competition system in the sense of Hirsch Il26ll27ll28ll 



(1.4) 



dvi{t) 
dt 



K 



Such systems has been intensively studied and many aspects are now well understood see for 
example ||IlB|3[l3[l5l[l3|2l|26l|2Zl|2i|29l|3^ and references therein. In particular, the 

existence of stationary solution and the asymptotic behaviour of the solution has been obtained 
in IIT2l|T3ll3Tll39] . Those systems are characterised by the existence of at least as many equilibrium 
states that the number of competitrng species (or genotypes) involved. In addition, the dynamics 
exhibit a competitive exclusion principle which state that the fittest species initially present will 
overcome all the other ones. 

When the mutations's matrix (/iy ) is non trivial, the system (|1.H does not fall into Hirsch's 
definition of competitive system and less results are known. If for reasonably smooth interaction 
functions, the existence of solution of (|1.H defined for all times is not an issue, the existence of 
a non trivial stationary solution and the analysis of the asymptotic behaviour are challenging 
questions. 
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Most of the known results concerns either particular interaction functions ^'^ for which the 
existence of steady states and their local stability are investigated l|4l [T5ll24l|29] or for some ODE's 
systems l|l.H where the mutation matrix /iy is considered as a small parameter. In the latter the 
system (|1.H is seen as a perturbation of (|1.4t and analysed using perturbative techniques ||6l|7l|8l. 
Recently, there has been an intense activity on continuous trait version of (|1.H where some of 
the techniques can be used to obtain the existence of locally stable steady states for l|l.Hl see for 
example l|3[Zl|8l|3|10lllllIl3|13l|l3|19l|3TJ|^ However, to our knowledge there is no 

results on global stability of the steady states for systems like (|l.H neither for its continuous trait 
version. 

In this direction, our first results concern the systems lll.H where the competition terms 
are independent of i. A tjrpical case is ^'^(w) = X^jLi ^'j which corresponds to a situation where 
inter- and rntra-species (genotypes) competitions are equals (i.e. blind and uniform competition 
between variants, see Lafforgue et al. Il32l ). In this situation one can show that there exists a 
unique positive stationary solution of (|1.H v e R^, which attracts all the trajectories initiated 
from any nonnegative and non zero initial data. Namely we show that 

Theorem 1.1. Assume that the interaction ^'^ is independent ofi and is satisfying the assumptions (|1.2t - 
(II. 3t , then there exists a unique positive stationary solution v to (|1.H . Moreover for any nonnegative 
initial datum w(0) not identically zero, the corresponding solution v{t) of converges to v. 

A case of particular interest is when the interactions i take the following form '^i{v) = 
12j'=i^j • This particular structure of interaction was initially introduced on a theoretical ground 
by Sole et al. ||44| to model the competition between viruses. Recently, this form of interactions 
has been used to explain experimental results of virus evolution llT6ll22ll37| . Sole et al. showed 
that Eigen's model of molecular quasi-species p^O] was to a large extent equivalent to the Lotka- 
Volterra competition equations under this assumption. 

For this type of interaction, we can rewrite the system as follows 



,1.5, if) 



1 ^ 



N 



For lll.SII besides the asymptotic behaviour of the solution obtained as an application of Theo- 
rem [TT| we can precise the speed of convergence to the equilibrium. Furthermore we can give 
an estimate of the time to reach near the equilibrium. In epidemiology, this type of informa- 
tions is of practical use for building tractable nested models. Nested models are a class of model 
which explicitly links the relationships between processes at different levels of biological orga- 
nization. They are often used to study the pathogen evolution by linking the disease dynamic 
of within- and between-host, see Mideo et al. Il36l . Their formalisation becomes more simple 
when the wi thin-host pathogen dynamic is faster than the between-hosts dynamic. Indeed, in 
such cases, using for example slow-fast reduction techniques commonly used in ecology 12), the 
within-host dynamic can be approximated by its equilibrium, see for example ||23l . 
More precisely, we show that 

Theorem 1.2. Assume that the interactions take the form 'i'i{v) = ^2^=1 '''j'^j then for any nonneg- 
ative initial datum v{Qi) not identically zero, the solution v{t) of | |1.5| | converges exponentially fast to its 
unique equilibrium. That is to say there exists two positive constants Ci and C2 so that 

\\v-v\\oc < Cie-<^^*. 

Next, we analyse the situation where the functions ^'i are not reduced to a single function. 
From a biological point of view, this situation could appears as, often, genotypes (species) exhibit 
particular association patterns (see for example ||45ll in the case of plant viruses). A way to model 
such phenomena is to take ^'^ of the form 'i>i{v) = X^jLi o^i^jTj'^j where a^j are crowing index 
||34ll25ll . This index is equal to 1 when the 2 species are distributed independently and ranges 
from (complete avoidance) to a large constant (near overlap) according to patterns of species 
association. 
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In this general context, our first result concerns the existence of a positive stationary solution 
for the system l|l.H assuming we have the following extra condition 

N 

(1.6) ^Mij<-^- 

This condition makes senses in our application framework since mutation rates (expressed as 
substitutions per nucleotide per cell infection) range from 10^'* to 10~^ for DNA viruses and from 
10~^ to 10"'^ for RNA viruses |43l whereas growth rates Vi are of the order of the imit or higher. 
For example, the overall growth rate of the RNA virus VSV was estimated to 0.6 virus /hour IflTl . 
Each cell infected by a single VSV particle produces from 50 to 8000 progeny virus particle/ cell 
infection ||47| . 

Under the extra assumption (II. 6t we prove 

Theorem 1.3. Assume that r.;, ^E'^ and fj,ij satisfy the assumptions (|1.2t -| [L3t . Assume further that il.6\ 
holds . Then there exists a positive stationary solution v to the system jl.W . 



Lastly, we obtain the convergence to the equilibria for some general interactions ^i. Namely, 
we show that 

Theorem 1.4. Assume that Vij'i/i and satisfy l|1.2|l - lll.3ll . Assume fyrther that M.6\ holds and 
\I/j(ii) = a(v) + eil^i{v) with Tpi £ Cf^^ uniformly bounded and a e Cl^^. satisfying il.2\ - fL^ is so 
thatW X € R^, Va(x) > 0. Then there exists eq so that, for all e < to, the positive stationary solution 
of HI. It attracts all the possible trajectories initiated from any non zero and nonnegative initial data. 

Note that Cl^^ perturbation of the particular interaction function a{v) = '^jLi '''j'^j satisfied 
the assumptions of the above Theorem. 

1.1. General remarks. 

Before going to the proofs, we want to make some general remarks and comments. 

Mathematical comments : First as remarked in [13|, we can interpret the system of equations dl.St 
as a discrete version of the continuous model 



(1.7) '-^-v[t,.) 



- 77 / r{y)v{t,dy) 



fj.{x, y)v{t, dy) - fj-{x)v{t, x) 



by posing for each sub population corresponding to a typical trait Xi G 



N 



v{t,x) = y2vi{t)S^., r{xi) = n, n{xi,Xj) = and = / n{x,y)dy. 

iR" 

Recently there have been a lot of works dealing with (|1.7|l and generalisation of it, see for 
example ||3|3[Zll|9l[10l[ni[l3[l3[lE[lS|311|33|l^ references therein. A large part 

of the analysis are concerning (|1.7|l in absence of mutation (ie yu = 0) or in the limit fi ^ 0. 
In the latter case, much have been done in developing a constrained Hamilton-Jacobi approach 
to analyse the long time behaviour of positive solution of this type of models see for instance 
||3l [l9l . Other approaches based on semigroup theory have also been developed to analyse the 
asymptotic behaviour and local stability of the stationary solution of l|1.7|l see [7'. ^ . Although 
some of the techniques developed in this two frameworks may be adapted to analyse the system 
l|l.H most of them fail when we try to prove the global stability of the stationary states. To 
tackle this difficulty we construct a set of Lyapimov fimctionals for the solution of l|l.H when 
the interactions 'fi are independent of i. Properly used these Lyapimov functionals enable us 
to analyse the asymptotic behaviour of the solution of jl-l} in this particular situation. From 
this analysis, we derive some consideration on the asymptotic behaviour of the solution of l|l.H 
in a general situation. The Lyapunov functionals are constructed in the spirit of the relative 
entropy introduced for linear parabolic operators in ||35ll . Even if our problem is nonlinear such 
type of relative entropy can still be constructed. These is an interesting new feature of nonlinear 
dissipative system. 
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It is worth noticing that a similar construction can be made for l|1.7|l giving us access to a simple 
way of analysing the asymptotic behaviour and the global stability of steady solution of il.7\ , see 

m. 

Along some of the proofs we notice that the existence of a steady state in Theorem 11.11 and 
Theorem II. 31 can be generalised to situation where the mutation matrix (/iy ) is not symmetric. 
Indeed, when the interactions vP^ are independent of i (Theorem ll.H , the construction of a unique 
stationary solution relies only on the Perron-Frobenius Theorem which holds true for non sym- 
metric matrices. However, for the general case (Theorem 11. 3t to obtain the existence of a steady 
state, we do require that the condition (II. 6t is replaced by 



(1.8) 



N , 



< 



Biological comments : First, we emphasize the biological interpretation of our result and particu- 
larly the role of the mutation term. Indeed, under biological compatible assumptions concerning 
the competition ipi, the mutation matrix (/iy ) and the growth rate r.i we have shown that, in 
sharp contrast with the classical results known for the Lotka-Volterra system without mutation, 
the mutation term deeply changes the dynamics of the system. On one hand, the mutation term 
stabilizes the dynamic of system by reducing the number of equilibrium up to a single equilib- 
rium and on the other hand, the mutation term precludes the competitive exclusion principle to 
occur. 

Second, we emphasize the biological relevance of relaxing some hypothesis regarding (i) the 
monotony of ^'i and (ii) the symmetry of the mutation matrix /iy to study viral demo-genetics 
d5mamics in structured hosts. To illustrate this point, let consider 2 patches pi and p2, and 2 virus 
genotypes vi and W2- The demo-genetics dynamics of the viral population in this system can be 
modelled by 



dw 
'dt 



Mw + Rw - ^'(w)w 



where w is the vector (vi^p^ , U2,pi , wi.pa : "2,p2 ) R, M and ^'(w) are the following matrices R := 



M := 



V 





1 - 









(-di 


d2 








fi fj. 










-di 
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^ di 





~d2) 
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In structured hosts some tissue types often act as virus "sources" and others are "sinks", creat- 
ing an asymmetry in the exchange between patches. In the above example this implies that di > 
d2 if the patch pi is a "source" and the patch p2 is a "sink", making the mutation matrix M non- 
symmetric. Moreover, in this example since the competition takes place only inside a given patch, 
we can check that the two monotone interaction functionals involved 4*1 (w) := ?'i^'i,pi + r2V2.p1 
and ^2{v) := rivi^p^ + r2V2,p2 do not satisfy the monotone properties l|1.2b . 

Such type of structure, inducing asymmetries in the exchanges and some weak interaction 
functionals, are expected to play a role in constraining or facilitating adaptive evolution of viruses 
in heterogeneous host environment ||38l . It is then relevant to study extension of our results in a 
context of more general assumptions on //^ and ^j. 
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Organization of the Paper : In Section |2] we start by establishing some preliminaries results about 
the system (|1.H . We derive the relative entropy identities in |3] Then in Section |4] we analyse in 
details the system (|1.H for a particular type of interactions and we prove the Theorem 11.11 We 
show Theorem 11.21 in Section |5] Finally, we are proving the Theorems 11.31 and 11.41 in Sections |6] 
andia 

2. Global facts on the system ijlTTIi 

In this section we establish some useful properties of solution of l|l.lt and prove the existence 
of a positive global in time solution of the system l|l.Hl , that for convenience we rewrite 

(2.1) ^ = A{^{t))v 

where A{'i'{t)) is the following matrix: 

/(''I - - Ml + Mil My 



A := 



\ Mij (rN - ^^Nit)) - ^iN + UnnJ 



with '^i{t) := and /i; = X^jLi Mij- Since the function "^i are locally Lipschitz the local ex- 

istence of a solution of (|2.1t is a straightforward application of the Cauchy-Lipschitz Theorem. To 
obtain a global solution starting with non-negative initial data, we need more a priori estimates 
on the solutions (wi)je{i,...,JV}- Let us first show that the system II2.H preserve the positivity. 

2.1. Positivity. 

Lemma 2.1 (Positivity). The system (|2.1b preserves the positivity. That is to says that ifv(0) is a non- 
negative initial value, then the solution v of \1A\ stays nonnegative. Moreover, for all i Vi{t) > Ofor all 
times te (0,T]. 

Proof: 

We argue by contradiction and assume that v changes its sign. Let ti e (0, T] be the first time 
where ~ for some i, v{t) > for all t e (0,ti) and there exists t > ti so that v{t) < 0. 

<i > is well defined since v changes its sign and Vi{0) > for all i, Vi{0) > for some i and 
linif^Q+ > 0. Now at the time ti, we have 

Therefore for all i Ui(ti) = = "^"'^^^^ and by the Cauchy-Lipschitz Theorem we deduce that for 
all i and all t G [ti, T] Vi{t) = 0. Thus u > in (0, T] which contradicts that v changes its sign. 

Next we show that Vi{t) > for all times t E (0,r]. From the above argumentation, v is 
nonnegative for all times. To show that v stays positive for all time, we can see that it is enough 
to show that 7V(t) = J2f=i W > for all times te [0, T] . Indeed, if there is <i € (0, T] so that h 
is the first time where for some i, Ui(ti) — then arguing as above we see that for all i and T > ti 
v,{t) =0 so Af{t) = Ofor all < > ti. 

Now, let us denote Q{0, 1) the following set 

Q(0,1) := |xeK^|^|x.|<l|. 
Since is locally Lipschitz, we have on Q(0, 1) 



i=l 



N 

(2.2) Wi,3Ki so that *i(x)<Kj^|xi 

1=1 
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Let Ko J2iLi ^ ^ positive constant so that 6 < min |l, We will prove that 

Af > S for all times in (0, T]. We argue by contradiction and assume there exists ti < t2 so that 
Af{ti) = 6 and Af{t) < 6 in (^1,^2]- Using that the Vi are non negative functions and l|2.2t , from 
our assumption we deduce that on [ti, and for all i, 

dv ^ 

>Vi {r,mn - KiAf) + ^ fJ'ljivj - Vi). 

By summing over all the possible i we end up with the following equation 

(2.3) ^>^f{r„^^n-^fKo), 

(2.4) M{h)=S. 

Using the comparison principle on the latter equation, we achieve on [ti ,t2], Af > Af where Af is 
the solution of the logistic type equation (12. 3L By construction we have Af > 6 for all times t > ti. 
Thus we get the contradiction 5 <Af <5. Hence TV > (5 for all t e (0, T]. 

□ 

Remark 2.2. Of the above proof, we also have a bound from below for Af. Namely, we have for 
alH 

AAW>min(l,^(0) ^" 



2.2. Existence of a global solution. 



2kq 



Next we show the following estimate an 

Lemma 2.3. Let v{t) he a solution of the Cauchy problem \2.1\ , then there exists two constants Co and 
Ci and a positive vector Vp so that v{t) < Coe^^^Vp 

Proof: 

From Lemma |2. II we know that the Vi are non negative functions. Therefore by \2.1\ we can 
see that the Vi satisfy 

(2.5) ^ < + ^ {vj - Vi). 

Let us denote M and R the two following matrices 

(-Ml + Mil \ 
R 
fJ-ij —^J■N + ^J'NN/ 

We can then rewrite the inequalities (12.51 as follows 

(2.6) ^<(i? + M)«. 

By choosing p, = sup^g j^y jii, we can see that i? + M + /2 is nonnegative matrix. Since R + M + 
jlld is also irreducible, by the Perron-Frobenius Theorem R + AI + fild posses a unique principal 
eigenpair {i^p, Vp) so that Vp is a positive vector, i.e. there exists {vp, Vp) so that Vp > and 

{R + M + jlld)vp = VpVp. 

Note that e^''''~f'>*\\v{0)\\^Vp satisfies 

du , „ , „ 
— = {R + M)u. 

So, by the Cauchy-Lipschitz Theorem, from | |2.6t we deduce that v < e'''p^'^'*||u(0)||oo^^p for all 
times t. 

□ 



fri 
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The existence of a global in time solution for the system (|2.1t is then a consequence of the 
Cauchy-Lipschitz Theorem and the above a priori estimates. 

Remark 2.4. By adapting some ideas in [39| and for a particular type of ^'i, we can derive an 
explicit formula for the solution of i\2.1} . Indeed when the interaction terms take the form ^'i (w) = 
Q^l'^) = SjLi '^i^i * then the solution Vi{t) can be expressed by the formula 

(g(i?+A/)t^(0))^ 



(2.7) v,{t) 



i + E;=i«.7o (e(«+*'Mo))^. 



To obtain this formula we argue as in {39l (Chapter 2 Section 2.1 ) and we start by introducing 
the fimctions 

TV 

a{t) ajVj and V^{t) := e^o 'i'v,{t). 

3 = 1 

We remark that the Vi satisfy the linear equation 



(2.8) ^Ml=r,V,+Y,ix,,[V,-V,). 



N 



dt 



Thus V^{t) := (e(^+^^)*w(0))^ and v^{t) is implicitly given by the formula 
Now let us evaluate the term -^o By differentiating e~ -^o one has 

, N 

_(e/„'«(«)<i«) ^ Q-(i)e/o ^ajyj(t). 

Therefore we have 
Hence 

(g(i?,+A/)t^(0))^ 



w,:(t) 



i + E;=i«.7o(e(«+^^Mo))^.d5 

□ 



3. Relative Entropy Identities. 



Here, we prove the following general principle which give us access to some useful identifies 
that we constantly use along this paper. 

Theorem 3.1. Assume that r^, {^lij), satisfies Hl.St . Let v and v be respectively a positive 

solution and a positive stationary solution of and let H hea smooth (at least C^) function. Then the 

function Hiv) J2tLi ^'fH (^j-^ satisfies 

N 

dt ' K 



,3.1) ^.-B(„) + li:.,H'f|)r,„, 
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where 



N 

Proof: 

By Hl.lt , for all i we have 



H { ^ 



Vi 



N 



(f) 













(3.2) 



Using that v is a stationary solution, we have for all ; 



1 ^ 



and we can rewrite the above equation as follows 

TV 



UVi ^ / Vj Vi \ 1 ^ 



By multiplying the above equality by ViH' I |^ J and by summing over all i we achieve 



1=1 

Thus we have 



_ 1 



dt K 



1 ^ /u\ ^ 



Vi I \ Vi V, 



Hence we have 



N 



since by symmetry of /i^ , 



N 



H{ ^] -H 

Vn J \ Vi 



□ 

As immediate corollary of the Theorem |3.1[ we have the following identities that we constantly 
used along this paper. 

Lemma 3.2. Assume that r^, (^y), satisfies (|1.2t - 1 |1.3| |. Lef v and v he respectively a positive solution 
and a positive stationary solution of l|l.H , then we have 

(i) 

d 
di 



/ N \ ^ N 

^ I " ;^ Xl*^*'*^^) " '^i{v))v,v, 

) 4=1 



(ii) 



Vi Vi \ I 



ij = l z=l 
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Proof: 

(i) and (ii) can be obtained straightforwardly from the Theorem |3.1| Indeed, by using H{s) — s 
in Theorem 13. li and by observing that from the symmetry of fiij, 

N 
*J=1 



0, 



we easily get 



dUjt) 
dt 



1 ^ 

— '^TiV^Vi, 



which proves (i) since for the function H{s) = s, we have 'H(t) = obtain (m) 

in a similar way by using the function H{s) = in the Theorem 13.11 and by observing that 

1 2 



N 
ij = l 



N 



4. The special case = «(.): 

In this section we analyse in details the asymptotic behaviour of a positive solution of l|l.H 
when the competition functional is independent of i, i.e for all i, '^i{v) = a{v) where a is a 
fimction from MJ^ K which satisfies l|1.2t - lll.3ll . As we expressed in Theorem ll.ll that we recall 
below, in this situation the system l|l.H has a unique positive stationary solution which attracts 
all the trajectories initiated from any nonnegative and non zero initial data. More precisely, we 
prove 

Theorem 4.1. Assume that for all i the function = a{.), then there exists a unique stationary 

solution V of jl.l^ . Moreover, for all nonnegative and non zero initial datum v{Q), the corresponding 
solution v{t) converges to v. 

To prove the Theorem, we first analyse the existence of stationary solution 

4.1. Study of the existence of equilibria. 

Recall that we look for a stationary solution of 



(4.1) 



dv . , , > > 
- = Aiait))v. 



Therefore if there exists a stationary equilibria for the system (|4.H the Vi must satisfies the 
following equations: 

(4.2) A{a)v = 

(4.3) a = aiy) 
where A (a) is the matrix 

A{a) := 

Note that we can rewrite the matrix A{a) := (i? — Id + M) where R and M are matrices 
defined in Section|2] 

Therefore, a solution v to l|4.2t is a solution to 

(4.4) ^M + R)v^(^^)v. 

Let us now establish some important property of the equilibrium. 

Lemma 4.2. Ifv is a nonnegative stationary solution of ^AAy then either v = Oorv > (i.e V z, > 0). 
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Proof: 

First, we observe that is a solution of the problem l|4.2t . Now let v be nonnegative stationary 
solution of (|4.H so that Vi = for some i. From (|4.4t we see that 

{M + Rv), = (-^) V, = 0. 

Therefore 

JV JV 

= TiVi + ^ ^^ij{vj - Vi) = ^ ^^ijVj. 

j=l 3 = 1 

Thus Vj = for all j since by assumption vj > and (/i^ ) is irreducible. Hence any nonnegative 
stationary solution is either positive or the zero solution. 

□ 

Observe that from the above Lemma and l|4.4ll , one can see that w is a positive eigenvector of 
the matrix M + R associated to the eigenvalue Now we are in position to prove that there 
exists a unique positive stationary solution to l|4.4t . 

Lemma 4.3. There exists a unique {a, v) solution of the equations (|4.2t and | |4.3t . Moreover v satisfies 
a{v) = fp. 

Proof: 

By choosing jl ~ supjgjj^ ^i, we can see that R + M + p. is nonnegative matrix. Since 
R + M + jlld is also irreducible, by the Perron-Frobenius Theorem R + M + Jlld posses a unique 
principal eigenpair [up, Vp) so that Vp is a positive vector, i.e. there exists {i^p, Vp) so that Vp > 
and 

(4.5) {R + M + flld)vp = VpVp. 

Moreover, the linear subspace associated to the eigenvalue Vp is one dimensional |46|. So without 
any loss of generality we can assume that X]i!=i("p)? = 1- From the equation Il4.5t we deduce 
that the vector Vp is a positive eigenvector of the matrix M + R associated with the eigenvalue 
Ap := {vp — fl). By construction one can see that Ap is the imique eigenvalue of the matrix AI + R 
associated with a positive eigenvector. A quick computation shows that Ap = (i/p — p.) > 0. 
Indeed, if not we have 

{R. + M)vp < 0. 

Thus for alH g {1, . . . , we have 

N 

n{vp)^ + ^n^j{{vp)j - (up)j) < 0. 

Let (wp)io := minig{i ..._^}(i)p)i then for {vp)^„ we have 

JV 

^lJ-tj{{Vp)j - {Vp)^o) > 0. 

i=i 

Since i? is a positive matrix we achieve the contradiction 

N 

Now from l|4.4t we deduce that there exists an unique positive a so that = Ap. Let us now 
construct our solution. Note that for any /i <E R, the vector piVp is also a solution to (|4.4b with the 
eigenvalue Ap. So to obtain a solution v to (|4.2t and (14. 3t we only have to adjust fi in such a way 
that a{fiVp) = a. This is always possible for a imique /i since q;(0) = 0, limp^.oo a(/^Wp) — +oo and 
a is an increasing fimction. 

□ 
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4.2. Convergence to the unique equilibrium. 

Let us look at the convergence of v{t) toward its equilibrium. Let us first establish some useful 
identities. 



Lemma 4.4. Let us denote {a, v) the stationary solution constructed above. Let v be a solution of (|4T 
then V satisfies the following identities 

(i) 



(ii) 



^ N N / \^ 2 ^ 



i—1 ^jj — 1 ^—1 

Proof: 

Since here for alH, ^'i = a, from Lemma [3.2l we deduce that 



E " ~ E M^j^^Wj ( ^ - ) + -^{a - ^(a(u) - a{v)) [v^f 



Thus (i) and (ii) hold true. 



□ 



From the above Lemma we can derive a useful Lyapunov functional. 

Lemma 4.5. Let v be the positive stationary solution of (|4.H . For any positive solution v of let us 
denote E, /3 and F the following quantities 

N N 

p{v) ■■^Y.''^^^ F{v):=\og 



£{v) 



Then F > log( g^ ^ - ) and for any positive solution v of l|4.1t we have 



N / s 2 

V, Vi 



d_ 

Proof: 

First let us show that F is bounded from below. By construction and using a standard convex- 
ity argument, we see that 0^{v) < (supj Vi)£{v). So from the monotonicity of the log, we conclude 
thati^(«)>log(^). 

Now let us show that F is non increasing. To do so let us compute J^F(v). From the definition 
of F we see that 

d_ dtSjv) ^ dtPjv) 

dt ^ ' £{v) ^ I3{v) ' 

where dt denotes Jj. By Lemmas I4.4l we have 

dt£{v) 1 ^ - - f'Vj v.iY 2 

dtf3{v) 1 ... 
^ = -(«-a(t)). 
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Thus, 



N , .2 

Vn Vi 



T,nv) = --1- V ^^,,v,v, ( ^ - ^ I < 0. 
dt E(v) f-^^ \v, vj 

□ 

Next we derive some a priori estimates on the solution v of (|4.1t from the previous Lemmas. 
Namely, we show that 

Lemma 4.6. Let vbea positive solution of (|4.1t , then there exists Ci so that 

£{v)+P{v) < Ci. 

Proof: 

Observe that from Lemma [4. 51 to obtain the bound it is sufficient to have a uniform bound on 
/3. Indeed, since F is decreasing in time, we have for all times 



Now recall that by Lemma l4!4l /3(v) satisfies the following equation 

Since a satisfies the assumptions l|L2t - (IL3t , there exists i?^, Cq, so that for all x G M^'^ \ 



(4.6) Cc,|^^xJ <a{x). 

Assume that for some t > 0, v{t) E R^'+ \ Qr^ (0) otherwise we are done since 

P{v) < { max Vi ] \vj\. 
\^e{l....,N} J fr{ 

Let S be 

S := {t e M+ |7V(t) > Ra}. 
From iA.6\ and using that the Vi are non negative, we see that for all t e S 

ka 



ait) > c^N{tf- > (^] , 



where v^^ax maxj^ji ^v} '^i- Therefore, on Zl we have 

Using the logistic character of the above equation, we can check that 

Piv{t)) < sup I /3(w(0)), max /3(a;), ( ^ 



We are now in position to prove the convergence of v toward its equilibrium. 



□ 



Lemma 4.7. Let {vi)i=i...N be the unique stationary solution of ji.li . Then for any non negative initial 
datum (wi(0))i=i...Ar not identically zero, the corresponding solution (wi(i))i=i...Ar of (|4.H converges to 
{vi)i=i...N CIS t goes to infinity. 



14 



JEROME COVILLE AND FREDERIC FABRE 



Proof: 

For simplicity we denote <, > the standard scalar product in M^. 

Now, since v and for all times t, v{t) = (wi, . . . ,vn) is a vector of M^, we can write 
v{t) := X{t)v + h{t) with for all t, < h{t), v >= 0. Substituting v by this decomposition in (|4.1|l , it 
follows that 

(4.7) X'{t)v + = Xit)A{a{t))v + A{a{t))h, 

at 

(4.8) ^ ^{a-a{t))X{t)v + A{a{t))h. 

K 

Therefore, we have 

< X'{t)v + >=< J^(^- + A{a{t))h, h > . 

Thus 

dh(t) , 1 d£(h) , , 

By following the computation developed for the proof of (ii) in Lemma [4.41 we see that 

dE(h) f hj hj 



J2 \^.-^.] a{m{h) 



z,j — 1 ^ ■' ^ 

Since > for all times, we will analyse separately two situations: Either £(/i(t)) > for all 
times i or there exists to £ K so that £{h{tQ)) = 0. In the latter case, from the above equation we 
see that we must have £{h{t)) = for all t > to and so for all t > to, we must have v{t) = X{t)v. 
Hence from | |4.8| | we are reduced to analyse the following ODE equation 

X'{t) = ^ia~a{Xm 

where 3 is the increasing locally Lipschitz fimction defined by a{s) a{sv). Note that since 
X{t) < v,v >= (3{v) > 0, we have X{t) > for all times t. The above ODE is of logistic t5^e with 
non negative initial datum therefore by a standard argumentation we see that X{t) converges to 
A > where A is the imique solution of 3(A) = a. By construction we have 3(1) = a, so we 
deduce that A = 1. Hence, rn this situation, v converges to v as time goes to rnfrnity. 
In the other situation, £{h{t)) > for all t and we claim that 

Claim 4.8. £ih{t)) -^Oast^ +oo. 

Assume the Claim holds true then we can conclude the proof by arguing as follows. From the 
decomposition v{t) = X{t)v + h{t), we can express the fimction (i{v{t)) by f3{v{t)) ~< v, v >= 
X{t) < v,v >. Therefore from Lemma l4!4l we deduce that 

(4.9) X'it) = ^(a - a(A(i)w + hit)))Xit). 

K 

Now by using £{h) 0, we deduce that /i as t ^ oo and from (14. 9t we are reduced to 
analyse the ODE 

X'{t) = ^(a - a{X{t)))X{t) + X{t){a{X{t)v) - a{X{t)v + h{t))) 
A 

= ^{a~a{X{t)))X{t)+X{t)o{l) 



where |o(l)| = \a{X{t)v) — a{X{t)v + h{t))\ < Cy/£{h) — > as t — > oo. As before we can conclude 
that X{t) 1 and v converges to v. 

□ 

Proof of ClaimlOl 
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(4.10) 



Since £{h{t)) > for all t, as in Lemma [4.5l we have 

£{h) 



^ 1 



1 ^ 



Thus the fiinction F := log (^pl^^-jy^ is a decreasing smooth function. 

First we observe that the claim is proved if there exists a sequence (t„)„gN converging to in- 
finity so that £{h{tn)) 0. Indeed, assume such sequence exists and let {sk)keN be a sequence 
converging to +oo. Then there exists fco and a subsequence {tn,^)keN of (t,i)„gN so that for all 
k > ko, we have Sfe > tn^ . Therefore from the monotonicity of F we have for all k > fco 



log 



ii3{v{sk))r 



< log 



By letting fc to infinity in the above inequality, we deduce that 

£{h{sk)) 



lim log 

k—¥oo 



= — oo, 



which implies that £{h{sk)) — > 0, since by Lemma \A£\ (d(v(th)))hcN is uniformly bounded. The 
sequence {sk)keN being chosen arbitrarily this implies that £{h{t)) — )- as t -> +oo. 

Let us now prove that such sequence (i„)neN exists. We argue by contradiction and assume 
that inf(g]{+ £{h{t)) > 0. Therefore from the monotonicity and the smoothness of F we deduce 
that there is cq e M so that 

F{h{t)) Co and ^F{h{t)) as +oo. 



Thus by Lemma [46] and | |4.10t it follows that 



(4.11) E^^.-.-^(|-|)'=0- 

From the a priori estimates of Lemma 14.61 there exists a sequence t„ ^ oo so that for all i 
hi{tn) — > hi. Passing to the limit along this sequence in the equation (|4.1H it yields 



N 



= ^ flijVjVi 



hj hi 



2 



By using the irreducibility assumption on the nonnegative matrix ^ij and the positivity of the 
quantities Vi, one can deduce from the above equality that we must have for all i and j 

hj hi 

Vj Vi 

Thus if we set A := we have h = Xv. So by using that < h,v >= for all time it follows that 
A = 0. Hence, we get the contradiction 



< ■mi£{h{t)) < £{h{tn)) 0. 



□ 
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5. A CASE OF INTEREST 

In this section we analyse more precisely the d3mamics of the solution v of l|l.Hl when the 
interactions '^i take the form {v) '^jLi ^j'^j ■ prove the Theorem 1 1 . 21 that we recall below. 

Theorem 5.1. Assume that the interactions take the form 'l'i(w) = ^2^=1 '''j'^j then for any nonneg- 
ative initial datum u(0) not identically zero, the solution v{t) of | |1.5| | converges exponentially fast to its 
unique equilibrium. That is to say there exists two positive constants Ci and C2 so that 

\\v-v\\oc.< Cie-^''. 

Before proving this Theorem we establish two auxiliary results that for convenience we present 
in two separate subsections. We prove Theorem 1 1.2 1 at the end of this section. 

5.1. Study of the evolution of the total population. 

Let us denote Af{t) = J2iLi total population. A straightforward computation shows that 

for the interactions ^'i of the form ^'^(f ) :~ ''"j'^j ^^e that J^{t) satisfies the equation : 

which written with the new variable P{t) = K — M{t) takes the form 

The d5mamic of the above equation is strongly related to the behaviour of a{t) and we have 

(5.3) |P(i)| = |Po|e^-t J■o"('')'*^ 

Lemma 5.2. For any nonnegative initial datum i;(0) not identically zero, M{t) converges exponentially 
fast toward its unique equilibrium K. Moreover M satisfies identically 

(5.4) UnTL{N,mn{t),Nrriax{t)} < N (t) < max{N,mn{t) , Nmax{t)} , 

where Mmin und Afmax are the solutions of the logistic equations: 



du , 



dt ^ V' 
u(0) = AA(0) 

with respectively ^~ = min{ri, . . . , rjv} and ^+ = niaxjri, . . . , r^} 
Proof: 

Assume for the moment that Il5.4t holds then the convergence exponentially fast to X is a 
straightforward consequence of (15. 3t . Indeed by (15. 4t we deduce that Rmin niin(A/^m, A/^iax) < 
a{s) < Rrnax ram{J\fmt„ Mniax) whcrc i?inax := niax{ri, . . . , ttv} and i?,„i„ := min{ri, . . . , r^}. 
Therefore a{s) > Co, since JVmin and Nmax converge to K. 

To obtain (|574t we investigate the following three cases, 7V(0) = K, N{Q) < K and 7V(0) > K. 

In the first case A/'(0) = K, we see that ( I5.4t holds true trivially, since Nmax = Mnin = A/'(t) = 
K for all t. Let us now investigate the two other situations. The argumentation in both situation 
being similar we expose only the case < K. 

In this situation, Mmin and Afmax being the solutions of logistic equations they are increasing 
functions. Moreover wa have Afmin{t) < Mmax{t) < K for all t. On another hand, by continuity 
oiN{t), there exists also tx > so that A/^(t) < K on [0, ti). Furthermore on (0, ti) we can see that 
M{t) satisfies the following differential inequalities 

(5.5) ^ > Rm^n^^ {l - ^ 

,5.6, f ^«".".^-(-^ 
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From the above differentials inequalities, by comparing M , Afmm and Nmax via the Cauchy 
Lipschitz Theorem, we obtain A/'min < A/" < A/^,nax for all < e [0,ti). Since A/'ma2;(ii) < -ftr,wecan 
bootstrap the above argument and show that (15. 4t holds true for all t. 

□ 

5.2. A useful functional inequality. 

Next we establish a useful functional inequality satisfied by vectors h E v-^ where denotes the 
linear subspace of M.^ orthogonal to v. 

Lemma 5.3. There exists Ci > so that for all h € v 



± 



Moreover Ci = A2 where A2 is the minimal eigenvalue strictly positive of the linear eigenvalue problem 



N - N 
V, 



Proof : 

Let T be the following Rayleigh quotient 

N 

Observe that the first part of the Lemma is proved if we show that 

(5.7) inf I{h) > 0, 

or equivalently 

(5.8) inf I{h) = inf I{h) > 0, 

h<£v^,e{h) = l hev^ 

since for all real /i, I{h) = I{p,h). 

To obtain I I5.7L we argue by contradiction and assume that inf^gg^ T(ft,) = 0. By | |5.8t we can 
take (/in)nGN a minimising sequence so that /i„ € v-^,£{hn) = 1 for all n. Since {x G |a; E 
v'^,£{x) = 1} is a closed boimded set, (/in)neN is uniformly bounded and we can extract a sub- 
sequence (/i„fc)feeN which converges to h e v-^ with f (/i) = 1. Passing to the limit along this 
subsequence, we see that l(h) = which combined with h G v-^ implies that h = 0. Thus we get 
the contradiction = £{h) = 1. Hence | |5.7| | holds true. 

Let Ci be defined by Ci = inf, £(,1)^1 let us try to compute Ci. By construction Ci 
is the result of a constrained minimisation problem. Therefore, by the standard optimization 
Theory 1,42 J , the minimizers must satisfy the following Euler-Lagrange equations 

(5.9) {D - M)h = Xh, 

(5.10) £{h) ^ 1, 

(5.11) <h,v>=0, 

where A e M is a Lagrange Multiplier to be determined and M and D are the following matrices 



M ■= 



D := 



n 1 V^-'V 

V ^Ej^ii^NjVjj 



By taking k > large enough, the irreducible matrix M — D + kid becomes non-negative 
and so M — Z? + kid is a symmetric positive definite matrix. As a consequence, the eigenvectors 
associated to the eigenvalue of M ~ D + kid form an a orthogonal basis of M.^ , i.e. = 
0^ El,. , where E,,. denotes the eigenspace associated to the eigenvalue Ui. 
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Since M — D + kid is non negative and irreducible by the Perron-Frobenius Theorem there 
exists a unique eigenvalue, says i^i associated with a positive eigenvector (f)i. Furthermore i^i 
is the largest eigenvalue and is algebraically simple. By a direct computation one can see that 

(M — D + kld)v = kv, thus vi ^ k and 0i = ^v. From the properties of vi we have Ey^ — lin{v), 
for aWi ^ l,Vi < vi — k and v-^ — E^^ . 

By construction, the A; — fc < are the eigenvalues oi M — D and we can see that for all 

h € v-^,£{h) = 1 we have 

I{h) = %{M - D)h > min{-AJ. 

Let A2 < be second largest eigenvalue of M — D and (f>2 an associated eigenvector By 
normalising (j)2 properly and since 4>2 & v-^ a straightforward computation shows that 

Ah) = *MD - M)(j,2 = -A2 - miii{-AJ. 
Hence, min^jg^i "^{h) = —^2- 

□ 

5.3. Asymptotic behaviour of the solution. 

We are now in position to obtain the exponential rate of convergence for any solution v of | |1.5t . 
Proof of Theorem ll.2l 
First we claim that 

N 

(5.12) Y,v, = K. 

1=1 

Indeed from Lemma [4.31 we deduce that v ~ — —Vp where Vp is the positive normalised 

eigenvector associated to the principal eigenvalue Xp of the matrix R + AI. A straightforward 
computation shows that 

N N N 

i—l i—1 i—1 



Since A„ > it follows that 



Thus 



. N 



Now recall that in the proof of Lemma [4.71 from the orthogonal decomposition v(t) ~ X{t)v 
h{t) we have 



By Lemma [5.3| we obtain 



"^^^^^ <-C,E{h) + ha-a{t))E{h). 



dt - ^ ' K 
Now arguing as in the proof of Lemma [4.5[ we end up with 



Therefore, thanks to Lemma [46l we deduce that 

' £{h{Q)\ 



(5.13) £{h) < I f/y^y, ) e-^i*/?2(«) < Cse-^^* 
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Finally we get the exponential rate of convergence, by observing that 

7V(t) V, 1 >= A(t) < u, 1 > + < /i, 1 >, 
which thanks to (15.311 implies that 

K </i,l> |AA(0) - ifle-i/X")'''* 

X(t) = ; \ ^ h — Z ■ 

< W, 1 > < V,l > <V,1> 

By Lemma [5.21 and the estimates | |5.13t , | |5.12t , ,using standard norm estimates, we have 

</i, 1> |A/'(0) -A'le-^./'o 



\V - U oo < 



< 



K K 



N8{h) 



□ 



6. The general case:The stationary solution 



In this section we investigate the existence of a positive stationary solution of \2.\\ under the 
additional condition lll.6t on the matrices R and M that we recall below, 

N 

This assumption has for consequence that the matrix R + M is positive definite. Indeed, we 
have 

N ^ N N , N 



-i— 1 "i^j — ^ ^tJ — 1 'i'-J — ^ 



> 

- 2 



Thus kern{R + M) ~ {0} and the matrix _R + M is invertible. Moreover from the last inequality 
we see that there exists positive constants co and Co so that 

(6.1) Co < M, u >j?,+j\/< £[u) < Co < U,U >B^+M, 

where < u,u >r+m:~ *u{R + M)u. 

Let S(ii) be the diagonal matrix defined by 

With this notation, a positive stationary solution of ^2.1} is then a non negative solution of the 
following problem: 

(6.2) (i? + M)v = E{v)v 

Note that when can be written as S(v) = a(v)Id, the construction of a positive solution 
has already been made in Section|4] So in the later, we will assume that cannot be written 
as S(w) = a{v)Id. It is worth mentioning that in this situation the method used in Section|4]does 
not work and we have to use another strategy. 

Let T be the following map 
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T 



V ^ Tv= {R + M)-'^[E{v)v] 

Since R + M is invertible, T is well defined and one can easily check that a positive solution of 
i6.2\ is a positive fixed point of the map T. To check that T has a positive fixed point we use a 
degree argument. 

Let ^'i(w)''' and E{v)^ defined by 

*f s^, + (1 ~ s)vl/i, {E^{v)),, := 

we consider the homotopy H e C([0, 1] x K^, M^) defined by 

i? : [0, 1] X ^ 

{s,v) ^ := + 

One can see that H{1, .) =T and H{0, .) ^ Tq where Tq corresponds to the map 

To : ^ 

Note that there exists a unique positive fixed point to Tq which can be constructed by arguing 
as in SectionlH 

The next step in this degree argument is to obtain for all s, a good a priori estimates on the 
fixed point of the map H{s,v), i.e. a good estimate on the positive solutions of the following 
problem: 

(6.3) (R + M)V = '^%V)V. 

In this direction we show the following: 

Lemma 6.1. Let V be a non negative solution of i6.3\ . Then either V = OorV>0 and there exists ci 
and Ci independent of s so that 

i 

Proof: 

To obtain that V is either positive or V = we can argue as in the proof of Lemma 12.11 So 
assume that there exists a io so that Vig = 0. By construction Vi„ is a minimum of the Vi and from 
the equation satisfied by Vi^ we get 

TV 

Therefore Vj = Vig for all j where i-iigj ^ 0. Since M is irreducible there exists j ^ iq so that 
Mioi 7^ 0- Let p := {k\Vk = 0} then iq and all j so that fii^j ^ belongs to p. In the previous 
argument, by replacing io by any k E p, we see that all j so that i^kj 7^ belongs to the set p. By 
iterating enough times the above argument and using the irreducibility of the matrix M we can 
see that p ~ {1, . . . , N} so Vi = V^p = for all i. Therefore a nonnegative solution of i6.3\ is either 
= or 1/ > 0. 

Now let us assume that V > 0. Recall that by iS.li there exists positive constants cq and Co so 
that for all ueR'^ 

Co <U,U >B+M< £{u) <Co<U,U >R+M ■ 

So for a solution V of < |6.3| | one has 



<V,V >R+M= V ^nV)V,^ >[s min + <i/i{V)) < V,V > 



CONVERGENCE TO THE EQUILIBRIUM IN A LOTKA-VOLTERRA ODE COMPETITION SYSTEM WITH MUTATIONS 21 



and we also get 

<V.V >R+M< Co (s max ^,{V) + *i(l^)) <V,V >r+m ■ 

Therefore we have 

(6.4) ^'i(F) < 

Co 

(6.5) i-<s max + l'i(y). 

Co ie{2....,N} 

Now thanks to the assumptions (|1.2t - (|1.3t made on the functions ^'i, there exists i?i , ci , fci and 
N positive constants Ki so that : 

/ N \ 

(6.6) Forall xeM^^+\gfl,(0), ci(^a;,j < *i(x), 

N 

(6.7) For all i, and for all x e Qr^^O), < \xj\. 
By combining I I6.4MI6.6L i6.5\ and (16. 7t we deduce that 



^ r 1 

A/" = y V, > mill <^ Ri, — 

^ I Co(ki + SUPj Ki 



□ 



6.1. Computation of the degree. 



We are now in position to prove the existence of a positive solution to the equation | |6.2t by 
means of the computation of the topological degree of T — irf on a well chosen set ft C 
Now we take two positive constants C2 and C2 so that C2 < ci and C2 > Ci where ci and Ci are 
the constants obtained in Lemma 16^1 Let be the following open set 



N ^ 

J2v.<cA 
1=1 J 



and let us compute deg{T — Id,Q,0). 

By Lemma l6.ll for all s e [0, 1] in H{s,v) — v ^ on dfl. Therefore using that H{., .) is an 
homotopy we conclude that deg{T - Id, Q, 0) = deg{H{\, .) - Id, f], 0) = deg{H{0, .) - Id, Q, 0). 
By construction, from Section l4?Tl one can check that deg{II{0, .) — Id, 17, 0) 7^ since the map Tq 
has an unique stable and positive fixed point. Thus deg{T — Id, ri, 0) 7^ which shows that T has 
a fixed point in Vl. 

□ 



7. The general case: Asymptotic Behaviour 

In this section we prove Theorem 11.41 That is to say under the extra assumption | |1.6t we 
analyse the asymptotic behaviour of the solution u(t) when for all i the interaction can be 
expressed like: '^i{v) = a{v) + ei/'i(w) with ipi uniformly bounded. To obtain the asymptotic 
behaviour in this case, we follow the strategy developed in Subsection 14.21 Namely, we start by 
showing some a priori estimates on the solution, then we analyse the convergence by means of 
a Lyapimov functional. For convenience we dedicate a subsection to each essential part of the 
proof. 
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7.1. A priori estimate. 

We start by establishing some useful differential inequalities. Namely we show that 

Lemma 7.1. Assume that ri,{^ij),'^i satisfies M3\ and (II. 6t . Assume further that = 

a{v) + eipi{v) with Tpi uniformly bounded. Then there exists eq > so that for allO < e < eq there exists 
a;+ G M^, w"*" positive and a positive real 7 so that 

(i) 



N 



i=i 

> ^(«(7'^+)-a(«))E 



(ii) 



dt 



< 



> 



N / 



(iii) 











E{v) 





< 



> 



1 -+-+ / ^'i 



+ _(a(7cI,+ )-a((i+)) 



Proof: 

First, we observe that (iii) can be straightforwardly obtained by combining (i) and (ii). So we 
deal only with (z) and (m). 

Let us denote a := e\\ip\\oo- Then since v is positive, from l|l.H it follows that 



dv. 
~dt 

dv- 



N 



< {r, + cr - a{v))vi + ^p-ijivj - f^), 

A 



Let u)~^ and be the stationary solutions of the corresponding equations 

dut 



dt 

du£ 
dt 



A 



(r; ~a-a{uj ))uj^ + ^ (w^- - ). 
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Now, let us fix e small enough so that ui^ exists. Then by arguing as in the proof of Lemma 
we obtain 

d ^ 1 ^ 

i=i ]=i 

d ^ 1 ^ 



Note that by Lemma l43l cj* are positive eigenvectors of the matrices R + M± aid. Thus ui are 
eigenvectors of the matrix R + A4 associated to the principal eigenvalue of R + AI . Since R + M is 
irreducible, the eigenspace associated to the principal eigenvalue is unidimensional. So, we have 
a;+ = 70;^ for some positive 7. Hence (i) and {ii) hold true since 

N / X^AT / \2 

Vi V, \ \ - _+_+ / Vj \ 



2 



□ 



Next, we derive some a priori estimates for the solution v for an interaction as in Theorem 
L4I We first prove some sharp a priori estimates for stationary solution of (|Llt . 



Lemma 7.2. Assume that ri, satisfies \1.2\ . (IL3II and M.6\ . Assume further that ^^(w) = 

a{v) + t'il)i{v) with Tpi uniformly hounded. Then there exists ci < Ci and ei so that for allO < e < ei 
and for any positive stationary solution of jl.l^ , we have 

N 
i=l 

Moreover, there exists tq so that for all e < ei , G Q^ro i^^)- Furthermore, for any nonnegative initial 
datum Vi{0) not identically zero there exists two constants C2(w(0)), (72(^(0)) so that for all < e < ei 
the solution satisfies for all t, 

N 

C2 < PMt)) ■.= J2ve,^Ve4t) < ^2- 
i=l 

Proof: 

Let us first observe that for e < eo by replacing by in (i) of Lemma [7. 11 we get 

N 



(7.1) 0<^ia{i^+)-aiv,))J2^iv,,^, 

1 ^ 

(7.2) o>^a{jti+)-aiv,))J2^tve,: 



From the proof of Lemma |431 we also see that 
(7.3) a{ij+) = K{Xp + cr), a{uj-)=K{\p~(j), 

where Ap is the positive principal eigenvalue of the matrix R+M. By using the monotonicity of a, 
we deduce that the maps ct i-)- are monotone. Moreover, we have < a(cj+) — a{ui^ ) < 2aK. 
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To obtain ci and Ci we argue as follows. Let us fix ci minjeo, — }■ Then by Il7.3t , for all 
e € [0, ei] we have 

(7.4) a{Co+)<a{u+)^^-^, a(c^-) >«(.;-) 



4 ' ^ 
Now thanks to the assumptions (|1.2t -< |1.3t satisfied by a, there exists Ra ,Ca,ka and Kq, so that 



/ N \ '^^ 

(7.5) Forall a; e R^^+ \Qii„(0), c^^x, <a{x), 

N 

(7.6) For all i, and for all x G Qr^ (0), a{x) < Ka ^ \xj\. 

So for e e [0, ei], by combining ||7lJ, lE3,|[73},lEl, EH} and we achieve 

(7.7) V,, < sup I ('^) ~ ,rA^:C,, 



(7.8) ^wej>min<jii'„,5^ 



i=i 

N 

. I „ I 

: ci. 



i=i 



To obtain a more precise estimate on v^, we argue as follows. Let us decompose AeW^ + 
where is orthogonal to Then by replacing by in (iw) of Lemma |7. 11 and using that 

— a{u!^)\ < 2aK we get 

N / \ ^ 



(7.10) a E I - ^ I 



2 



From ^7.7} and using the functional inequality, Lemma [5.3[ we deduce that 
(7.11) SiK) < i^. 

Combining the latter estimate with ^7.7} and the positivity of aj+ and v^, we have the estimate 



C'l Na 
A. < — 1 + 2/ 



where uj is the stationary solution with e = 0. Let Rq :~ Ci ^1 + , and choose ei smaller 

if necessary to have i?o < 26*1. Next consider the set (5i?„ (0) c Q2Ci(0)- From the above estimates 
on and h^, since a is Lipschitz continuous in Q2C1 (0) there exists kq so that for all e. 



|a(AeW+) - a{XeU)+ + h,)\ < nnyjN£{h,), 



which combine with (|7.1H enforces 

(7.12) \a{Kuj+) - a{KCj+ + K)\ < 2CiKQy 

Thus from ([7l), and l(7l2ll we deduce that 



(7.13) |a(cS+) - «(A,w+)| < 2Ka + 2CikJ^^ < C^/^, 
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with C independent of e. Observe that by construction there exists two positives constants < 
to < 1 < ti so that for all e < t\, there exists G (toj'^i) so that aj+ = l^Qi. Recall now that 
by assumption Va > 0, then the real map s n- ai^suj) is smooth (C^(M)) and increasing. It 
is therefore an homeomorphism in R+ and a local diffeomorphism in R+. So by the Inverse 
Function Theorem, we deduce that for all s, t e (0, C), 

k 

\s -t\< —\a{su!+) - a{tid+)\, 



(-0 



where 



C:= sup{-, ^^ I 1 + 2/ 



k := 



1 

iin,se(o,c)(< V(q!(sw)),w >)' 



In particular, we have 
which combined with i7.13\ enforces 



k 

|s - 1| < — |a(sw+) - 



|Ae-l| < -CV^. 



Hence from the latter estimate and l|7.1H we have for all e < ei, f;^ G (5p^(<D+). 

Next, we derive an uniform upper boimd for f3{ve) when e £ [0, ei]. In the sequel of this proof, 
for convenience we drop the subscript e on ij. 

First, we observe that by Lemma [3.2| we have 

\i=i / 1=1 

1 \ ^ e ^ 

— - a{t)) ) ^ v^Vi + -jr'^iMv) - i^^{v))v^Vi. 



' 1=1 i=l 



Since the fimctions V'i are imiformly boimded, we achieve 



d ( ^ \ ( \ \ ^ 



ViVi 



By using dZU and ||Z!D, it follows 



d I v-^ _ \ 1 7AAp , , v-^ 



Again using that a satisfies the assumptions l|1.2|l - ||1.3l l, there exists Ra,Ca, ka so that for all x e 

M^'+\Qfl,„(0), 

(7.14) c„(^x,j <a{x). 

Now, let us assume that for some t > 0, v{t) e M^'+ \ (5_r„ (0) otherwise the proof is done since 
we have 

Biv) < I max Vi | > \vA < CiRa, 
V'e{i J fr{ 

where Ci is the boimd obtained above. Let E be the following set 

T,:={teR+\M{t)> Ra}. 
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From l|7.14|l , l|7.7|l and using that the are non negative, we see that for all t e S 

■ 1 



Therefore, with cq Ca (^) we have for t e S 



dt ~ K \ 4 

Using the logistic character of the above equation, we can check that 



I3{v{t)) <snpip{v{0)), sup I3{x), 
Thus, by using Il7.7t , we achieve for all e e [0, ei] 



4cn 



P{vc{t)) < C2 := sup |/3(«(0)), (7ii?„, C, (^^^ 

To obtain the lower bound for /3(we) we argue as follows. First let us observe that by taking ei 
smaller if necessary, since ui'^ — > a) as e — )• and e (5(^^(<D+), there exists a positive constant 
Co independent of e e [0, ei] so that for all stationary solution we have 

min Vei > Cq. 
ie{i,...,N} ' 

Now by Il7.8t , for all e e [0, ei] we have 

N N 

(7.15) (^(Ve) > min > Cq w, 

ie{i,....N} j—' j—; 



e,2 ■ 

1=1 i=l 



Therefore to obtain an uniform lower bound for it is enough to obtain an uniform lower 
bound for := ^c,*- From l|l.Hl by summing over all i and by using the definition of 5* and 

the boundedness of the tpi we deduce that A/^ satisfies the following inequality 

Thanks to (|7.H and H7.4t , we have 

By reproducing the argumentation of the proof of Lemma [2.1l and by using Remark l2.2l we can 
check that 

(7.16) ^.„)>„,„{i,^,m}, 

where ki denotes the Lipschitz constant of the function a in the unit cube. Hence, by collecting 
(I7.15t - Il7.16t we achieve for all e < ei and all t > 0, 

□ 

Remark 7.3. Note that from the above argumentation, using the Logistic character of the equa- 
tions, we can get that for all t <t\ and all initial data w(0) > 0, there exists to so that for all t > to 
we have 
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Lastly, we obtain some uniform control on a continuous set of homeomorphisms 

N 



where u e U c M^'+. 
Namely, we show that 

Lemma 7.4. Assume that ri, {fiij),'i'i satisfies jl.li , il.3\i and (11. 6t . Assume further that ^i(u) = 
a{v) + eipi{v) with a € C/^^ satisfying (|L2t , jl.Si and V'i € Cl^^ uniformly bounded. Then there exists 
€2 and To > so that for for all e < £2 and for all stationary solution of l|l.H we have 

C3 < <C3< 2(73 < +To). 

Moreover there exists and fc > so that for all e < e^we have for all stationary solution of (|LH and 

t,SG(0,l+To) 

Proof: 

Recall that from the proof of Lemma [7!2] for all e < ei, for any stationary solution of (|1.H w^, 
we have a(cj^) < a(we) < a{oj^). So, for all e < ei and for all stationary solution of (|1.H we have 

where a = e||-0||oo- Let us fix 62 := minjei, ^^jj'^y — }. From (|7.4t and Lemma [7.2l we get 

.-,:^2^<i.,(l)<52|^^:C3, 

for aU e < £2 and for any stationary solution of (|1.H . 

By using to the monotonicity of the map a ^ ui^ and Lemma [7.21 we can choose £2 smaller if 
necessary to achieve for any £ < £2 and any stationary solution v^, 

•so'^e^ < soVe < sqw^ for all sq > 0. 

The latter inequalities imply that we have for all £ < £2 and for any stationary solution v^, 

$B,(soWe) > a(sow7i) - 

Let us fix So such that 

This is always possible since a is monotone increasing and lim^_).oo ce{iiu!~_^ ) — +00. By construc- 
tion. So > 1, since ^' r,^ is monotone increasing and we have 

$c,(so) > 2$o.(l). 

Let us denote tq := sq ^ 1- 

Now since for each i the function satisfies l|1.2[l , il.3\ and that a and the ipi are C;q^(K^), 
the function ^'d, (s) := X^i^i is C;^q^(R) and monotone increasing. Therefore, for any 

stationary solution of (|1.H v^, the fimction vj/^^ is a M+ homeomorphism. 

Next, we check that for a fixed e and a fixed stationary solution Uj, the homeomorphism is 
a diffeomorphism on (0, 1 + to) — > ^'c^ ((0, 1 + tq)). Thanks to the Inverse Fimction Theorem, 
to show that ^Pu^ is a local diffeomorphism it is sufficient to prove that for all s € (0, 1 + tq), 
^'v^ (s) ^ 0. By a straightforward computation we have: 

N 

(*sj'(s) =< Va(sue),ue > +£^ < V^^{sv^),v > vli. 

1=1 
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For all s e (0, 1 + tq), and for any stationary solution £ Q(i+ro)Ci (0) we have: 

\y^p^isv,)\< sup |VV'»(w)| := Ci, 

''6Q(i+xo)<3i(0) 

ya{sVf ) < sup Vq!(u) := C2 > 0, 

■"eQ(i+xo)Ci(0) 

Va{sv^) > inf Va(w) Ca > 0. 

■"eQ(i+.o)Ci(0) 

Therefore, by Lemma [7.21 we deduce that 



(*.j'(5)>^C3-eCiC?. 

By choosing e < £3 := min |e2, 2JVCiC^ ^^'^ ^' ^ (0' 1 + '''0) and for any stationary 

solution of HI. Hi Ve, 

From the the latter a priori bounds, we see that for all e < e^, s E (0,1 + tq) and for any stationary 
solution of (|l.lt Ve, we get the following a priori estimate 



C3 



Hence, we deduce that for all e < £3, s, t e (0, 1 + tq) and for any stationary solution of 
?;e, we have 

|s-i| < fc|$s^(s)-$s.WI, 

with 

, _ 2N 

□ 

7.2. Asymptotic Behaviour. 

We are now in position to obtain the asymptotic behaviour of the solution (t) as t goes to 
+00 for e € [0, £*], where e* is to be determined later on. 

Let us fix e € [0, ei] where ei is obtained in Lemma 17.21 and let be a stationary solution of 
(|1.H . For simplicity we denote <, > the standard scalar product in R^. 

As in the proof of Lemma [4.71 we start by observing that since ^ we can write v^{t) := 
A(t)wc + h{t) with for all t, < h, >= 0. In the sequel of this subsection, for more clarity in the 
presentation we drop the subscript e on v and v. 

First, we note that from this decomposition we can derive the following equalities: 

N 

(7.17) X <v,v >^ y^^ViVj, 

(7.18) ^ = ^'-2AA' 

at at 

By Lemma Oand l(7l7t , we can check that for all e < ei and for all i > 0, we have 

(7.19) § < Mt) < 
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Similarly, by using Vi, we deduce from 

£{v) = £{h) + \^£{v) that 

(7.20) E{h)<^[N + 1), 

^1 

(7.21) j:M<^{i+^ 

, ci \ Cl 

By (l7.19t , ll7.20t , l|7.2lV Lemma [7.2l and by using the Lipschitz regularity of and the Cauchy- 
Schwartz inequality we can check that for some constant C > independent of e and v 

N 

(7.22) - < c^jm, 

i=l 
N 

(7.23) - < C^£{h). 

i=l 

Indeed by I l7.19t - ll7.21ll and Lemma [7l2l Xv and h are uniformly bounded and we have 

N N / N \ 

i=i 1=1 ) 

where Ki are the Lipschitz constant of t/i^ on the set (5/j„ (0) with i?o := {l + . From 

the decomposition of v and by using the Cauchy-Schwartz inequality, we get 

N I N N \ N 

(7.24) K'^^^^) - ^'(^w^ + < s ii - m\ E + E I'^^i E I'^'^i^^ 

1=1 \ j=l 3 = 1 I i=l 

(7.25) < NK^Sih)£{v) [l 1 - A(t) I ^£(5) + ^/£(}ij 

where R := supjgj]^ ^j. Ki. Therefore, by ||7.20|| and Lemma [7.21 the inequality Il7.23ll holds true 
for some positive constant C independent of e < ei. A similar argumentation holds to get l|7.22t . 
Next, we show that 

Lemma 7.5. There exists e* < miiijeo, ei, £3}, so that for all e < e*, E{h^{t)) ast +00. 



Assume the lemma holds true, then we can conclude the proof of Theorem IL4l by arguing as 
follows. 

By combining ^7.171 and Lemma [3. 2 [ we achieve 



d ^ 1 " 

(7.26) \'{t) < V, w >= — E = ^ E 



dt^ K 

2 = 1 1=1 



Now by using £{h) 0, we deduce that h ^ {) ast ^ 00 and from the latter equality we are 
reduced to analyse the ODE 



i=l 

where 

N 

0(1) :=E(*'(^W^)-**W)"^"''' 
1=1 

which by (|7!22l and |[7!23l satisfies 

C(l + e) 



N 



(7.27) |o(l)| < ^ ,7 '^ V^W^O as t 
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By Lemma l7!2l for all e < e* < ei we have C2 < J2iLi ^i'^i — ^2 with C2 and C2 independent of 
e. Therefore by l(7l7t ^ < X(t) < ^ for all t > 0. Thus A satisfies an ODE of the form 

which can be rewritten 

(7.28) A'W = -^(0 + 0(1) -$(AW)), 

where := and $(.s) By construction, * e Cl^c(^) and is 

monotone increasing. 

The above ODE is of logistic type with a perturbation o(l) with a non negative initial 
datum. Therefore, by a standard argumentation, we see that X{t) converges to A > where A is 
the unique solution of ^'(A) = 6. By construction, we have ^(1) = 0, so we deduce that A = 1. 
Hence, X{t) — > 1 and we can conclude that v converges to v. 

□ 

Let us now turn our attention to the proof of the Lemma [7.5 1 
Proof of Lemma [731 

First, let us denote := — ^'i(u)- Then, by combining l|7.17t , | |7.18t and Lemma [3. 2 [ we 

achieve 

^ N ^ N 

(7.29) X'{t) <V,V>^—J2 ^r^r 



and 



dt^ K 

i = l 2 = 1 



i,j = l ^ ■' "' ^ 2=1 4=1 



i,j=l ^ ^ ^ ^ i=l 

Therefore using the definition of ^'i and with the notation 7^ := tpi{v) — ipi{v), we have 



d£[h) ^ _ ^ _ ^\ ^ _ ^^^^^^^^^ 



dt ^ V w, w, / -ftT 



2e^ 26A(t)^ 

;^2^72^; + — ^2^72«2^2, 

1=1 i=l 



which implies that 



N 



(7.30) ^<-f^^,^v,v^{h_h^j +:^[a{v) + 2eU\\oo-a{v))£{h) 



2eXit)^ 

y ^jtVihj. 



K 

i=l 



By M.29\ , using the definition of '^i we also have 

d ^ 1 ^ e ^ 

Z— 1 2—1 i — 1 

1 ^ 

- ;^('^(^) ~ ^ellV-IU - a{v)) X ViVi. 
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By Lemma I7l2l for all e < ei, C2 < I3{v) — J2iLi ^i'^^i < ^2 with C2 and C2 independent of e. So we 
have 

I log(/?(^;(^))) > ].{a{v) - 2e||ViU " 



which combined with | |7.30t leads to 



d£{h) 
dt 



< 



N 

E 



hj hj 



dt 



log{l3'{vm£{h) 



K 



2e\{t) 
K 



N 



By using the functional inequality. Lemma 15751 and rearranging the terms in the above inequality 
we get 



(7.31) 



dSjh) 
dt 



K 



N 



£\h) H ^— 2_^-/iVihi, 



where Ci (e) is the second largest eigenvalue of some associated linear problem. 

Now, we estimate the last term of the above inequality. Recall that by | |7.25l l we have 



N 



1=1 

By combining the above estimate and <I7.3H , we achieve 
(7.32) ^ - £ih)^^ log(/?2(«(t))) < + eCs) £ih) + eC,\l - Xit)\^m, 



where Cf 



and Ca 



2C2NR 

K ■ 



Recall that Ci (w) is an eigenvalue, therefore for e small enough, say e < £4, since v G Qc^i'^^) 



and w+ — > a; as e one has Ci (w) > 



C'(m + ) 



With the latter estimate and by choosing e < min(ei , £4) smaller if necessary, we get 



(7.33) ^ - £{h)j^ \og{P\v{i))) < ^^^l^£ih) + eC,\l A(OI 

The proof now will follow several steps: 

Step One: Since by |[7l9t we have |1 - A(t)| < 



NC2 



for all t, we claim that 



Claim 7.6. For all e < niin(ei, £4), there exists to so that for all t > <o ^t's have 



^/£{h)<2e[^ 

C2 



NC2 



Indeed, by H7.17|l , and Lemma [7!2l we have 
(7.34) 



'-fl - £ih)l logi^Mm < + eC4 



NC2 



and we can check that there exists to > so that ^j£{h{tQ)) < 2- 



Ci(i.+ ) 



. If not, then for all 



t > y/£{h{t)) > 2 c[{ci,+')' dividing | |734t by y/£{h) and rearranging the terms, we 

get that 



(7.35) V^^log 



<-^i^V^+£C4 



NC2 



NC2 
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Thus F(t) :— log (;a4^7(t)y) is a decreasing function which by Lemma l7.2l is boiinded from below. 



Moreover, we have for alH > ^j£{h{t)) > 2- 



e4C4 


1+^' 




^1 . 



Therefore F converges as t tends to 



and ^ 0. Thus for t large enough, we get the contradiction 



< V£{h)-log 



£{h) 



dt V/3^(«W) 



NC2 









=1 . 



Let E be the set S ;= |t > to\\/S{h{t)) > 2 — c,(cd+)^ j- Assume that S is non empty 
otherwise the claim is proved since ^ > 1. Let us denote t* := inf E. So at t* we have 



ciCi 1 



NC2 



yJS{h{t*)) = 2 Ci(ii'+)^ • Again, by dividing| |7.34t by y/£{h) and rearranging the terms, on 

the set S we have 



< 0. 



Thus log 



WW)) 



is a decreasing fimction of t for all i G E. By arguing on each connected 



component of E, we can check that for all t G E 



1 



C2 CiiCj+) 

Therefore by using Lemma [7.2l f or t > t* we have 

"1 _|_ iYgi 



e4C4 



C2 Ci{iv+) 



Hence, since ^ > 1 we get for all t > to, 



1 + 



NC2 



C2 Ci(cS+) 



□ 



Sfep Two: First, we define some constant quantities: 



(5o :=2 



1 



7VC2 



C2 Ci(cS + ) ' 

, NC2' 



(7.37) 
(7.38) 

where the constants C, k are respectively defined in l|7.22t ,| [7!23t and in the Lemma l7!4l 
By the previous step, we have that for all e < e* , all t > to 



C3 



We claim that 

Claim 7.7. For all e < e* there exists tfS„ such that for all t > t^Sg 
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Proof: 

First, we can check that for e < e* there exists t* so that for all t>t* 

\1-X{t)\ < 2kCeSo. 
Indeed, by Il7.28ll and l|7.27t , we have 

(7.40) A'(t) = -^(0 + o(l)-$(AW)), 
with 0(1) < Ce5o for t > to- Let X±ceSo be the solution of the ODE 

(7.41) X'^a.soit) = ^^|y^(*«^(l) ± ^'^0 - *«.(A±c..o W))- 
Since e < 4^^, for t > to we have 

f^<$,.(l)±C.<5„<^. 

Therefore A±cei5o ~^ X±ceSo where X±ceSo are the unique positive solutions of ^'o, (A±cc(5o) = 
$5,(l)±Ce5o. 

Thanks to the strict monotonicity of ^5^, we also have 

X-2CeSo < X-CeSo < XceSo < X2CeSo7 

where X±2CeSo are the unique positive solutions of 4'^^ (A±2Ce<5o ) = (1) ^ 2Ce5o- 
Since e < e* < ^^^^^ , it follows that 

I < *,,(!) ±2C6<So<^ 
and therefore, by Lemma l|7.4t , we have 

(7.42) < A±2Ce5o < 1 + T-o. 
Now, recall that for t > to, X{t) satisfies 

A'(t) < + Ce6o *..(A(t))), 

X'it) > ^^L($,^(i) _ Ce6o - ^,Am))- 
Ktyv) 

Thus, by a standard argumentation, we can show that for t > t* we have 

X-2CtSa < X{t) < X+2CeSo- 

Therefore, for t>t* we have 

|1 - X{t)\ < SUp{|l - A_2Ce5o|, |1 - X+2CeSo\}- 

By (|7.42t and Lemma IZil since e < e* < £3, we deduce that for t>t*, 

|1 - A(i)| < fcsup{|*c,(l) - ^vA~^-2C.So)\, - *C.(A+2Ce5„)|} < 2fcCe(5o. 

From the latter estimate, by using < |7.33t , we see that for t > t* 

(7.43) ^ - £ih)^^ log(/?2(i>(t))) < -^^l^£{h) + 2e'C,kCdoVW)- 

at at 4 

By following the argumentation of Step one, we can show that there exists t^So such that for 
t > teSo we have 

2kCe^S^ 



d 

Hence, we have for all t > t^Sa 
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since eSo < j^. 

□ 

Step Three: Since for all t > t^Soi 

by arguing as in the proof of Claim [771 we show that for all e < e* there exists t s^so that for all 

t > t^So 

By reproducing inductively the above argumentation, for all e < e* we can construct a se- 
quence (t„),igN so that for all t > t„ we have 

Hence, for all e < e* we deduce that 

lim S{h{t)) 0. 

□ 
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